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ABSTRACT 



Aims. We present a new list of thorium and argon emission lines in the visible obtained by analyzing high-resolution (R = 1 10, 000) 
spectra of a ThAr hollow cathode lamp. The aim of this new line list is to allow significant improvements in the quality of wavelength 
calibration for medium- to high-resolution astronomical spectrographs. 

Methods. We use a series of ThAr lamp exposures obtained with the HARPS instrument (High Accuracy Radial-velocity Planet 
Searcher) to detect previously unknown lines, perform a systematic search for blended lines and correct individual wavelengths by 
determining the systematic offset of each line relative to the average wavelength solution. 

Results. We give updated wavelengths for more than 8400 lines over the spectral range 3785-6915 A. The typical internal uncertainty 
on the line positions is estimated to be ~10 m s -1 (3.3 parts in 10 s or 0.18 mA), which is a factor of 2-10 better than the widely used 
Los Alamos Atlas of the Thorium Spectrum (Palmer & Engleman 1983). The absolute accuracy of the global wavelength scale is 
the same as in the Los Alamos Atlas. Using this new line list on HARPS ThAr spectra, we are able to obtain a global wavelength 
calibration which is precise at the 20 cm s level (6.7 parts in 10 10 or 0.0037 mA). 

Conclusions. Several research fields in astronomy requiring high-precision wavelength calibration in the visible (e.g. radial velocity 
planet searches, variability of fundamental constants) should benefit from using the new line list. 

Key words. Atomic data - Techniques: spectroscopic - Atlases - Catalogs 



1. Introduction 

The wavelength calibration of astronomical spectrographs in the 
visible is usually performed using arc lamp spectra such as those 
produced in hollow cathode lamps. Thorium is today the most 
often used element due to its numerous spectral features over 
the whole visible and near-infrared domains. Other advantages 
of thorium include its mono-isotopic nature and the absence of 
hyperfine structure, which lead to narrow, highly symmetric line 
profiles. For all these reasons, ThNe and ThAr hollow cathode 
lamps have become the standard for wavelength calibration of 
astronomical spectrographs. 

The widely used reference for accurate Th wavelengths is 
the Atlas of the Thorium Spectrum obtained with the McMath- 
Pierce Fourier Transform Spe ctrometer (FTS) of the Na tional 
Solar Observatory at Kitt Peak dPalmer & Engle man! 1983L here- 
after PE83). The FTS scans from which this atlas was derived 
have a spectral resolution of about 600,000. The atlas gives the 
position of about 1 1,500 lines between 3000 and 1 1,000 A. The 
quoted wavenumber accuracy ranges from 0.001 crrT 1 for the 
strongest lines to 0.005 cm -1 for the faintest ones. In velocity 
units at 5500 A, this corresponds to 16-82 m s _1 (5.3-27 parts 
in 10 8 or 0.29-1.50 mA). 

While the wavelength accuracy of the PE83 atlas will be high 
enough for most astronomical applications, it starts to be a lim- 
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iting factor when aiming at very high-precision radial velocity 
(RV) measurements. Over the past decade, the need for such 
measurements has considerably increased following the discov- 
eries of the first extrasolar planets using Doppler velocimetry. 
The improvements of this technique over the years have led to 
the detection of planets having the mass of Neptune and be- 
low. This require s a Doppler precision of 1 m s _I or less (e.g. 
lLovis et aill2006l) . To carry out such observations successfully, 
all aspects influencing the RV measurement, from the properties 
of the source to the instrumental design and the data reduction, 
have to be carefully taken into account. As far as we are con- 
cerned, we have been particularly involved in the development 
and operation of the HARPS instrument, a cross-dispersed, fiber- 
fed echelle spectrograph mounted on the Europe an Southern 
Observatory 3.6m-telescope at La Silla, Chile dMavor et al.l 
120031) . HARPS is operated under vacuum in an isolated, strictly 
temperature-controlled environment and is therefore extremely 
stable and well suited to high-precision spectroscopic work. We 
have investigated many aspects regarding instrument calibration 
and data reduction, and identified several domains where im- 
provements could be made. One of these is the wavelength cali- 
bration process. 

Two main techniques are commonly used today to achieve 
high-precision RV measurements: the so-calle d iodine cell 
(Butler et al.l 1 19961) and simultaneous ThAr dBaranne et al.l 
1996) techniques. While the former uses a superimposed iodine 
spectrum as absolute wavelength reference, the latter relies on 
ThAr calibration spectra (this is the one we use on HARPS). In 
this case, the stability and precision of the wavelength solutions 
derived from ThAr spectra is a crucial aspect of the calibration 
process. It turns out that the precision of the Th wavelengths as 
given in PE83 is no longer adequate for the wavelength calibra- 
tion of high-resolution spectrographs if a RV precision of 1 m s _1 
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is to be reached. Indeed, errors in the PE83 atlas, which are due 
to random noise in the FTS scans, translate into systematic off- 
sets of the lines with respect to the average wavelength solu- 
tion in the spectral chunk to be calibrated. The residuals around 
the wavelength solution are no longer dominated by photon shot 
noise, but by these systematic offsets which can be one order of 
magnitude larger than random noise, as for example in the case 
of HARPS. As a result, wavelength solutions will not be well 
constrained and will show unstable behaviour against adding or 
removing Th lines in the fitting process. 

The limited precision of the PE83 wavelengths led us to the 
idea of creating a new atlas of Th and Ar lines based on HARPS 
ThAr spectra. Contrary to Fourier transform spectrometers, grat- 
ing instruments have complicated dispersion relations and their 
overall wavelength scale cannot be determined from the mea- 
surement of a single well-known transition. However, they are 
usually much more sensitive, which should make it possible to 
use such instruments to reduce the random noise on the PE83 
wavelengths while keeping the absolute wavelength scale given 
by the FTS. HARPS spectra are indeed of sufficient quality to 
permit significant improvements in the relative precision on line 
positions. Moreover, these spectra will allow us to identify thou- 
sands of unknown (mostly faint) lines, which are not present in 
the PE83 atlas due to its lower sensitivity. Finally, this would 
give us the opportunity to perform a careful treatment of blended 
lines, which represent a major problem when aiming at high- 
precision wavelength calibration. 

In this paper we describe the new line list we obtained us- 
ing HARPS ThAr spectra. The different steps of the process are 
detailed in Sect. 2, and the results presented in Sect. 3. A dis- 
cussion about the potential of the new atlas, its advantages and 
limitations, follows in Sect. 4. 



2. Experimental method 

2.1. Instrumental setup 

The light source for the Th spectrum is a standard ThAr hollow 
cathode lampQ operated at a current of 9.0 mA. The lamp win- 
dow, made of UV glass, has a diameter of 37 mm. The pressure 
of the fill gas is 800 Pa. In such lamps, a discharge is established 
between the anode and the Th cathode through the carrier gas 
(Ar in our case). Collisions between the accelerated ions and the 
cathode sputter off Th atoms from the cathode. Collisions in the 
plasma in turn excite Th atoms, which emit light by decaying to 
lower states. As a result, a rich spectrum of Th and Ar emission 
lines is produced. 

The beam from the ThAr lamp is injected into an optical 
fiber carrying the light to the spectrograph, which is located in- 
side a vacuum vessel in an isolated room below the telescope. 
Light is dispersed on a large R4 echelle grating and the echelle 
orders are cross-dispersed by a grism. The entrance slit (i.e. the 
fiber output) is then re-imaged onto a 4k4k CCD mosaic, made 
of two 2k4k EEV CCDs with a pixel size of 15 fim. The spectral 
format consists of 72 echelle orders covering the whole visible 
spectral range from 3785 to 6915 A, except the interval 5306- 
5339 A which falls within the gap between both CCDs. The 
spectral resolution is R — 110,000 and the average dispersion 
0.015 ApixeL 1 (820 m s pixeL 1 in velocity units). This leads 
to a sampling of ~3.3 pixels per resolution element. 
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1 The lamp was manufactured by Cathodeon/Heraeus (type 
3UAX/Th, serial number B 13773). 



Fig. 1. HARPS instrumental stability and precision on the mea- 
surement of line positions. The CCD position of a strong Th line 
is measured on 64 different ThAr lamp exposures taken over a 
period of one month. The measured scatter shows that the spec- 
trograph did not drift during this period within the error bar of 
less than 0.002 pixel. 



2.2. Characteristics of HARPS ThAr spectra 

In the following we will be using a series of ThAr lamp expo- 
sures obtained in December 2004 as part of the standard calibra- 
tion plan. The exposure time was 40 s for each frame. ThAr spec- 
tra have been extracted and flat-fielded using the standard echelle 
reduction procedures implemented in the HARPS data reduction 
software. The basic data products we are starting from are the ex- 
tracted two-dimensional (pixel vs. order) spectra. These contain 
thousands of Th lines, which are unresolved at the HARPS reso- 
lution and intrinsically highly symmetric. Moreover, the instru- 
mental profile itself is very close to Gaussian. In the following 
we therefore use simple Gaussian functions to fit the position of 
the lines on the extracted spectra. 

Line intensities extend over the whole dynamic range of the 
CCD. For the strongest unsaturated lines, the photon noise on 
the line position is about 0.002 pixel, corresponding to 1.6 m s _1 
on average (5.3 parts in 10 9 or 0.029 mA). On the other hand, 
the position of the faintest detectable lines can be measured with 
a precision of about 0.2 pixel (160 m s _1 ). However, in the case 
of very strong lines, photon noise is actually not the limiting 
factor. We estimate that CCD inhomogeneities, such as pixel 
size variations (up to a few percent), inhomogeneous intra-pixel 
sensitivity or CCD block stitching errors, will limit the preci- 
sion at which the true position of a line can be measured to 
about 0.01 pixel, corresponding to 8 m s _1 (2.7 parts in 10 8 or 
0.15 mA). 

To minimize photon noise, the most straightforward option 
would be to combine many ThAr spectra and measure averaged 
line positions. However, this is only feasible if the instrument 
remains sufficiently stable over the time duration of the ThAr 
sequence. This condition is usually not fulfilled under normal 
terrestrial conditions because variations in atmospheric pressure, 
temperature changes and mechanical flexures induce instrumen- 
tal drifts of tens to hundreds m s _1 over timescales of minutes to 
hours. In the case of HARPS, instrumental stability is such that 
the spectrum does not drift by more than 1ms" 1 over timescales 
of days to months (see Fig.[T). This allowed us to combine a se- 
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Fig. 2. Closest-neighbour separation between lines in the final 
list. 



ries of 64 ThAr spectra taken over one month to measure line po- 
sitions, thereby considerably reducing photon noise. For the vast 
majority of the lines, the measurement uncertainty could actu- 
ally be reduced to about 10 m s _1 (3.3 parts in 10 8 or 0.18 mA). 
A quick comparison with the numbers quoted for the PE83 atlas 
in Sect. 1 shows that line positions can be determined 2-10 times 
more precisely using HARPS ThAr spectra. The use of the PE83 
atlas leads to a dispersion of the residuals around wavelength so- 
lutions of 50-70 m s _1 due to the errors in the input wavelengths. 
The creation of a new line list would allow us to bring this num- 
ber down to about 10 m s _1 . This is the main motivation for the 
new atlas. 

2.3. Line-fitting process 

At the beginning of the process, we developed an algorithm to 
systematically search for lines in the ThAr spectra in order to in- 
clude several thousands of lines which are not listed in the PE83 
atlas because it is less sensitive than our HARPS spectra (except 
in the bluest part). The use of a much larger number of lines, 
even if most of them are quite faint, has a stabilizing effect on 
the wavelength solutions, especially in the red part of the spec- 
trum where strong Th lines are scarce. Thanks to the newly de- 
tected lines we were able to increase the total number of lines 
by -50%. 

This has a major impact on the occurrence of blends. ThAr 
spectra are actually so rich that truly isolated lines are very rare, 
even at a resolution of 1 10,000. Blends represent a major prob- 
lem for high-precision wavelength calibration because they in- 
duce shifts in the centroids of the components which are of- 
ten much larger than photon noise. Fig.[2]shows a histogram of 
the closest-neighbour separation between lines. The histogram 
peaks at ~4.5 km s _1 , which corresponds to ~1.5 times the spec- 
tral resolution. This obviously means that there remain many 
undetected blends within the spectral resolution. However, we 
automatically rejected lines that appear single in our spectra, 



but actually contain other components resolved by PE83 (at 
R = 600, 000). In this way we eliminated most of the danger- 
ous blends. The effects of unidentified, faint blending compo- 
nents within the HARPS resolution are less severe since they 
will mostly induce shifts comparable to, or smaller than, photon 
noise. To deal with blended lines partially resolved by HARPS, 
we developed an algorithm to simultaneously fit neighbouring 
lines with multiple Gaussians. The algorithm determines which 
lines cannot be fitted independently and groups them with their 
closest neighbours to make a multiple fit. This allowed us to keep 
in our list blended lines whose components are only ~ 1 resolu- 
tion element apart (2.7 km s~'), provided the components are of 
similar intensity. 

While the occurrence of too close blends was the main cause 
for rejecting lines, about 100-200 strongly saturated lines and 
those in their neighbourhood had to be removed as well because 
they were impossible to fit properly. At the other extremity of 
the intensity scale, some PE83 lines were too faint to be fitted, 
especially in the bluest part of our spectra. 

After completion of the list of lines that could be properly 
fitted, we ran the line-fitting algorithm on each of the 64 ThAr 
spectra of our series. We then computed for each line its mean 
position, rejecting outlying measurements using sigma-clipping. 
Lines which showed an unusual number of rejected measure- 
ments were removed from the list for being unstable to the 
fit. Reasons for unstable behaviour include blended components 
with large intensity differences, faint lines at the detection limit, 
pollution by neighbouring strongly saturated lines, possible in- 
trinsic instability of a few contaminant lines, etc. The resulting 
list is our starting point for the wavelength calibration process. 

2.4. Determination of updated wavelengths 

The basic principle for the determination of updated wave- 
lengths is extremely simple: a wavelength solution is fitted 
through the measured line positions (in pixels) using their PE83 
wavelengths. The residuals to the fit will not be dominated by 
random noise from HARPS spectra, but by systematics from the 
PE83 atlas. Corrections to the PE83 wavelengths are then sim- 
ply given by the residuals to the wavelength solution. Obviously, 
wavelength solutions themselves have to be sufficiently well 
constrained and match the true optical dispersion relation as 
closely as possible. However, the systematics we want to cor- 
rect were originally nothing but random noise in the FTS scans 
from which the PE83 atlas was derived. They are therefore ran- 
domly distributed as a function of wavelength and do not have 
any large-scale impact. The best wavelength solution will then 
still be the one which minimizes^ 2 , i.e. the sum of the (properly 
weighted) squared residuals around the fit. 

We emphasize that this method will yield updated wave- 
lengths which are internally more precise than in PE83, but the 
global wavelength scale itself will be the same as in PE83, since 
we made no attempt to improve the absolute wavelength accu- 
racy. However, this has no impact on high-precision RV mea- 
surements for example, since these are essentially differential 
measurements. 

Our ThAr spectra are split into 72 spectral orders that overlap 
each other to some extent. This means that a significant fraction 
of Th lines are detected twice (on two consecutive orders). If we 
now fit an independent wavelength solution on each order, we 
will find two different systematics for a given line in the overlap 
region. This is particularly unsatisfactory because it is precisely 
towards order edges that wavelength solutions are most poorly 
constrained. We therefore chose to couple the wavelength solu- 
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tions of all orders by imposing one and the same systematic for 
all lines appearing twice. In practice, we first had to determine 
which lines are detected twice by assigning a preliminary wave- 
length to all of them. This was done by fitting independent wave- 
length solutions on all orders. We then classified lines as be- 
ing unique or duplicate. Concerning wavelength solutions them- 
selves, we found that third-order polynomials A(x) = 2^ =0 akX k 
give an adequate analytical representation in the sense that the 
residuals around the fit are randomly distributed as a function 
of wavelength (no large-scale trends). Moreover, supplementary 
degrees of freedom would often increase the reduced^ 2 and are 
therefore not desirable. 

The coupling between all orders generates a large linear 
least-squares problem where lines contribute different types of 
equations depending on their classification: 

- Duplicate PE83 lines: two equations of type 

2] a k x) - AAi = Ay>E83 (1) 

- Unique PE83 lines: one equation of type 

2] a k x) = /t,-,PE83 (2) 

- Duplicate unidentified or Ar lines: two equations of type 
£a*x?-^ = (3) 

- Unique unidentified or Ar lines: no equation 

In the above equations, the at are the unknown polynomial 
coefficients of the wavelength solution for the corresponding or- 
der, Xj is the position of line ; in pixels, A/l, is the unknown sys- 
tematic offset of line i, /1,',pe83 is the PE83 wavelength of line i 
and Ai is the unknown wavelength of line i (for unidentified or 
Ar lines). We consider Ar lines as having unknown wavelengths 
because they do not appear in PE83 (which used a ThNe source) 
and because they are sensitive to physical conditions in the lamp 
(see Sect. 3). The individual systematics AAi are unknowns to be 
determined, but they must have a global property that we can use 
to introduce more constraints on them in the least-squares min- 
imization: they should be normally distributed around zero. We 
can therefore add a further type of equation (one for each AAi): 

AAi = (4) 

This has to be understood in a least-squares sense: these 
equations will ensure that, in the global solution, the properly 
weighted (see below) A/i, will have a mean value close to zero 
and a variance close to one. 

We now need to assign the correct weight to each equation 
so that we can solve our weighted linear least-squares problem. 
Uncertainties for each equation will be computed differently de- 
pending on equation type. For types Q] and [3] uncertainties are 
simply given by random noise from HARPS spectra, derived 
from the combination of the 64 individual spectra. They are com- 
puted as the quadratic sum of photon noise and detector-related 
noise (0.01 pixel). For most lines these error bars could be re- 
duced below 20 m s _1 . For equation type [2] uncertainties are 
given by the quadratic sum of HARPS random noise and the 
PE83 wavelength uncertainty for the line under consideration 
(see below). Finally, equations of type [4] will be weighted with 
the PE83 uncertainty alone. 

To determine the PE83 wavelength uncertainty for each line, 
we first computed the residuals of the PE83 lines relative to 
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Fig. 3. Relation between line intensities and the dispersion of 
residuals around wavelength solutions when using the PE83 at- 
las. This gives the errors on the PE83 wavelengths as a function 
of the measured line intensity on HARPS spectra. In velocity 
units at 5500 A, errors range from 20 m s _1 to ~160 m s _1 (6.7 
to 53 parts in 10 8 or 0.37 to 2.9 mA). 



the wavelength solutions independently obtained on each spec- 
tral order. We then plotted these residuals as a function of the 
line peak intensities, which were previously corrected for the 
wavelength-dependent HARPS instrumental efficiency. The re- 
sulting plot is shown in Fig. [3] Residuals are given in cm -1 be- 
cause the PE83 uncertainties expressed in these units are ex- 
pected to be almost independent on wavelength (see PE83). 
There is a clear relation between line intensities and residu- 
als, which actually just reflects the photon noise on the FTS 
scans from which the PE83 atlas was derived. We fitted a power 
law to these data to obtain an analytical expression relating the 
PE83 uncertainty to the measured line intensity on HARPS spec- 
tra. This relation is in very good agreement with the error bars 
quoted in PE83 (see Sect. 1). 

The weighted linear least-squares problem can now be 
solved. We have in total 7648 equations and 2082 unknowns. 
The ^-minimization yields best-fit values for the coefficients of 
the wavelength solutions and the systematics AAi, from which 
updated wavelengths can be obtained for all PE83, unidentified 
and Ar lines. The reduced^- 2 is equal to 1.57. Several factors 
can explain its somewhat high value, among which a slight un- 
derestimation of the detector-related noise on HARPS spectra, a 
non-perfect analytical model to describe wavelength solutions, 
and the presence of undetected blends that introduce small shifts 
in line centroids. We checked that varying slightly the error bars 
has a negligible impact on the final results. Consequently we 
preferred to keep our best error estimates rather than artificially 
increase them. 
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New vacuum wavelength 
[A] 


Uncertainty 
[A] 


Wavenumber 
[cm- 1 ] 


Intensity 
[e-] 


Previous wavelength 
[A] 


Identification 


3785.650454 


0.000149 


26415.53974 


7273 


3785.649743 


NoID 


3786.675812 


0.000129 


26408.38692 


11987 


3786.675314 


Thll 


3787.458325 


0.000120 


26402.93078 


24284 


3787.457347 


Aril 


3788.714137 


0.000160 


26394.17923 


3050 


3788.714472 


ThI 


3789.009459 


0.000170 


26392.12203 


2584 


3789.010941 


ThI 


3789.435335 


0.000157 


26389.15595 


2707 


3789.435169 


Thll 


3790.185751 


0.000527 


26383.93118 


6700 


- 


7 


3790.243720 


0.000114 


26383.52765 


126158 


3790.243943 


ThI 


3791.160081 


0.000128 


26377.15049 


4572 


3791.158240 


NoID 


3791.343954 


0.000146 


26375.87125 


2491 


3791.344838 


ThI 


3791.431959 


0.000119 


26375.25903 


16543 


3791.432308 


NoID 


3791.496302 


0.001561 


26374.81143 


1711 


- 


7 


3791.871484 


0.000113 


26372.20181 


68571 


3791.871327 


ThI 


3792.296894 


0.000140 


26369.24344 


2497 


3792.297101 


ThI 


3792.373700 


0.000149 


26368.70939 


2076 


3792.373871 


Thll 


3792.594893 


0.000126 


26367.17151 


4127 


3792.594391 


ThI 


3793.807442 


0.000114 


26358.74422 


13282 


3793.806899 


ThI 


3794.893018 


0.000118 


26351.20398 


4508 


3794.893246 


ThI 


3795.228217 


0.000125 


26348.87661 


2699 


3795.228593 


Thll 


3795.775429 


0.000112 


26345.07807 


11036 


3795.775842 


ThI 



3. The new line list 

The new list gives updated wavelengths for 8442 lines between 
3785 and 6915 A. The data, available in electronic form at the 
CDS, contain 6 columns, as shown in Table Q] The first col- 
umn gives the new vacuum wavelength of the lines as obtained 
in this work. The second column contains the estimated inter- 
nal uncertainties on line positions. The third column gives the 
wavenumbers corresponding to the new wavelengths. Peak in- 
tensities, corrected for instrumental response, are listed in the 
fourth column. The fifth column gives the previously published 
wavelengths when available. For Th lines, this is the PE83 
wavelength. F or Ar lines, wavele n gths have b een taken from 
INorlenl (U973I) and IWhaling et all (Il995l 120021) . For unidenti- 
fied lines, the value was set to Finally, the last column 
gives line identifications when available. 'ThI', 'Thll', 'ThIII', 
'NoID' and 'Contam' are all identifications from PE83. The first 
three are self-explanatory. 'NoID' stands for lines that could not 
be identified as known Th transitions by PE83, but are prob- 
ably due to Th. 'Contam' stands for iden tified contamin ants. 
'ArF and 'ArlF identifications are from INorlenl (119731) and 
Wha ling et al.l d!995l [2002). Previously unidentified lines are 
designated with '?'. It should be noted that the line intensities 
given in the table are only indicative since flux ratios between 
lines may exhibit variations from one ThAr lamp to the other. 
In particular, the Ar-to-Th flux ratio is very sensitive to physical 
conditions in the lamp. 

We made some comparisons between several lamps of the 
same manufacturer to check whether wavelengths themselves 
might be affected by changing physical conditions in the lamps. 
Although a detailed analysis is beyond the scope of this paper, 
two main conclusions can already be drawn. First, Th lines turn 
out to be individually stable in wavelength at least within the er- 
ror bars quoted in this paper. Second, Ar lines exhibit wavelength 
shifts from lamp to lamp amounting to a few tens of m s _1 . This 
behaviour is probably the result of pressure shifts, to which th e 
Ar atom is known to be quite sensitive (see Whal ing et al.l l2002). 
As a consequence, the Ar wavelengths given in the atlas have to 
be taken with caution. In general, we strongly recommend not to 
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Uncertainty on line position [m s~'] 

Fig. 4. Distribution of uncertainties on line positions in the new 
atlas (shaded histogram) and the PE83 atlas (plain histogram). 
Note the logarithmic scale on the vertical axis. 

use Ar lines to perform wavelength calibration if high precision 
is needed. 



4. Discussion and conclusion 

4. 1 . Comparison with previous line lists 

If we merge the best availab l e Th and Ar line lists in the literature 
dPalmer & Englemanlll983t INorlenl 1 19731 IWhaling et aDll995l 
|2002|) . the total number of lines in the interval 3785-6915 A 
amounts to 7321. However, an important fraction of these can- 
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Fig. 5. Difference in line positions between the new atlas and 
the PE83 atlas as a function of wavelength. Only lines with a 
PE83 error below 50 m s are shown. Both atlases have the 
same absolute wavelength scale. 

not be used when working at a resolution of 1 10,000 because of 
blending of individual components. Actually, the number of us- 
able lines from these lists amounts to about 3900. The new atlas 
contains ~4000 previously unidentified lines and ~500 known 
lines which could not be used before because of blends with un- 
known components. In total, the new list more than doubles the 
number of lines available for wavelength calibration. 

As far as precision is concerned, we have briefly described 
the error distributions in PE83 and the new atlas in Sect. 2. To 
illustrate this more completely, Fig. H] shows both distributions 
on the same plot. The gain in precision is clearly seen and fully 
justifies the creation of the new atlas. 

We mentioned in Sect. 2 that the new list has, by construc- 
tion, the same absolute, global wavelength scale as PE83. This 
is illustrated in Fig. [5] which shows the difference in line po- 
sitions between both atlases. As expected, no large-scale trend 
can be seen. The mean offset between both lists amounts to only 
0.08 m s" 1 (2.7 parts in 10 10 or 0.0015 mA). The scatter itself 
mainly reflects the errors in PE83. 

4.2. Precision of wavelength calibration 

As a final step of our work, we implemented a new wavelength 
calibration algorithm in the HARPS data reduction pipeline in 
order to use the new line list. We now routinely obtain resid- 
uals around the fit which have a weighted rms dispersion of 
about 10 m s _I (3.3 parts in 10 s or 0.18 mA), as expected from 
the combination of photon noise and uncertainties in the input 
line list. This means that the wavelength solution is precise at 
a level of a few m s _1 locally, and at the unprecedented level 
of 20 cm s _1 (6.7 parts in 10 10 or 0.0037 mA) globally (i.e. for 
example when computing a radial velocity using the whole spec- 
tral range). This can be demonstrated in the following way: we 
choose one observation of a star and reduce it many times using 
different wavelength calibrations, i.e. different ThAr exposures. 
This will yield each time a slightly different radial velocity for 
the star. The dispersion of these RVs is a measure of the pre- 
cision we have on the wavelength solutions. Fig. [6] shows the 
results of this test. Using a series of 30 ThAr exposures, we ob- 
tain a raw dispersion of 24 cm s _1 . Correcting for instrumental 



Fig. 6. Test of the HARPS wavelength calibration accuracy. 
The radial velocity of a star is computed from a single stel- 
lar exposure using 30 different ThAr exposures as calibration 
frames. The obtained RV dispersion is a measure of the preci- 
sion achieved by the wavelength calibration process. 



drifts at the level of ~10 cm s _1 occurring during the series, this 
number reduces to 20 cm s _1 . 

As a conclusion, we find that the quality of the new line 
list permits significant improvements in the wavelength calibra- 
tion process in general. A number of research fields in astron- 
omy requiring high-precision wavelength calibration will cer- 
tainly benefit from using the new atlas. Apart from radial veloc- 
ity searches for extrasolar planets, we can also me ntion the work 
on the variability of fundamental constants (e .g. iMurphv et al.l 
120031: iLevshakov et al .1120051: IChand et al.l2006l) . where an accu- 
rate wavelength scale is needed to properly measure line shifts 
possibly caused by varying physical constants over cosmological 
timescales. 

We recommend using the new list for calibrating all spec- 
trographs working at a resolution R < 110,000. Obviously, the 
lower the resolution, the more adaptations will be needed to deal 
with blended lines. For each blend, the expected shift of the main 
component can be estimated using the positions and intensities 
of all components given in the input list. This value can then 
be compared to the random noise on the main component (or 
to some tolerance value) to decide whether or not the line can 
be safely used for wavelength calibration. An example of such a 
procedure illustrati ng how to select suita ble lines at lower resolu- 
tion is presented in lMurphv et al.l (l2007al) . At resolutions higher 
than 1 10,000, the choice of the input line list depends on what is 
exactly needed. The new atlas should be sufficient to guarantee a 
precise wavelength calibration provided the number of detected 
lines is high enough. If more lines are needed, then it is neces- 
sary to include lines from PE83 that were removed from the new 
atlas because of blends. In this case, great care should be taken in 
selecting lines to make sure they are not blended with unknown 
components. 

4.3. Future developments 

Although ThAr lamps have several qualities as wavelength cal- 
ibrators, they also have a number of drawbacks. At moderate 
spectral resolution, blends represent a major problem which may 
prevent high-precision wavelength calibration. Furthermore, the 
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high dynamic range in line intensities makes it sometimes dif- 
ficult to find a good compromise between, on the one hand, the 
minimization of strongly saturated lines and, on the other hand, 
the availability of a large number of lines. 

Finally, the finite lifetime of lamps and possible aging ef- 
fects affecting Th wavelengths at the ~1 m s _I level make 
these lamps unsuitable for future instruments aiming at cm s _1 
precision. This is for example the case for CODEX (COsmic 
Dynamics Experiment), a project for a high-resolution spec- 
trograph t o be installed on the European Extrem ely Large 
Telescope (Pasquini et al. 2006; Grazia n et al. I I2007I) . New cali- 
bration systems have to be developed to ensure a wavelength ac- 
curacy and stability at the 10~ n level. This has to be compared 
to the accuracy level of a few parts in 10 8 achieved by PE83. 
The needed "ideal" calibrator would produce a dense pattern 
of equally spaced lines with precisely known and stable wave- 
lengths. Moreover, the lines should all have about the same in- 
tensity. All these requirements are obviously challenging, but ef- 
forts are presently being made t o adapt the laser frequ ency comb 
technology to these needs (see Murp hvet al.ll200 7b). The pos- 
sibility to couple the laser comb to an atomic clock would per- 
mit to reach the required very high stability level. There is little 
doubt that, once fully adapted, such a calibration system would 
find applications in a number of different domains. 
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